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Abstract 



All real solutions of the Lane-Emden equation for n — 5 are obtained 
in terms of Jacobian and Weierstrass elliptic functions. A new family of 
solutions is found. It is expressed by remarkably simple formulae involving 
Jacobian elliptic functions only. The general properties and discrete scaling 
symmetries of these new solutions are discussed. We also comment on their 
possible applications. 



is a classical equation of mathematical physics. Introduced in 1870 by Lane Q], 
and studied in 1907 by Emden 0, it was originally used to describe the mass den- 
sity distribution inside a (spherical) polytropic star in hydrostatic equilibrium (in 
that case, the Lane-Emden equation follows directly from the Poisson equation for 
the gravitational potential and the assumption of hydrostatics). Shortly after Em- 
den's work was published, Plummer [3] proposed a model of the distribution of 
stars in globular clusters that was based on solutions of Eq. ([I]). Since that time, 
the Lane-Emden equation and its generalisations were applied in many different 
branches of physics, ranging from astrophysics to kinetic theory and quantum me- 
chanics (see, e.g., |4] and references therein). 

There are 2 well known solutions of the n - 5 Lane-Emden equation, i.e., 



1 Introduction 



The Lane-Emden equation 




(1) 




(2) 



The first of them, 



1 



(3) 



m = ± 



i 



is due to Schuster and Emden O. The second, 

m = ±A=, ( 4 ) 

which is singular at £ = 0, belongs to a larger class of solutions valid for all Lane- 
Emden equations with n > 3, that is 



0(0 



2(n - 3) 



Because the Lane-Emden equation (|2]) is invariant under the transformation 
S:6>(£) -> 0(^M)/VI, /i 6 K+, Schuster's integral Q yields, in fact, a whole 
one-parameter family of solutions. Solution Q is a fixed point of 5 . 

It was believed for many years that finding other solutions to Eq. Q 'is com- 
plicated and involves elliptic integrals' (Chandrasekhar |5]). Then, in 1962 Srivas- 
tava 10 found another solution that can be written in a compact form, namely 

sin (In V?) 

m = ± , - — , (5) 



^-24 sin 2 (in V?) 



where again the scaling transformation gives a whole family of solutions. This 
solution is also singular at £ = 0, but it can be used in composite stellar models 0. 
It was rediscovered in 1977 by Sharma 00, who searched for non-elliptic solutions 
of©. 

About a decade ago Goenner and Havas [4] found the following expression for 
another class of solutions 

m = Cl (6) 



+ ^(ln(^)/2; 4/3, -8/27 + 16^/3)' 



where p denotes the Weierstrass elliptic function, and B and c\ are integration 
constants (cf. Eq. (30) of 0). This simple result was somehow overlooked by 
researchers — see, e.g., a recent book by Horedt [9], where the aforementioned 
opinion of Chandrasekhar is repeated. 

In this paper we derive all remaining real solutions of Eq. ([2]) (new solutions are 



given by Eq. ( 10 1 in Sec. 2.1). We also show that for a certain range of integration 
constants, solution (|6]) can be obtained in terms of Jacobian elliptic functions. Be- 
cause of many different physical applications of Lane-Emden equations, we will 
not discuss any particular boundary conditions here. 



2 Solutions 

Equation Q can be transformed to an autonomous form by the following substitu- 
tion: 6 = zl an d f = — ln£. The resulting equation is 



d 2 z 1 / 4 n 



2 



and its standard integration yields 

■>2 



where C denotes an integration constant (see e.g. ll5l). 

Further analysis depends on the factorisation of the polynomial w(z) = -z 6 + 
3z 2 + C appearing in Eq. (T7j). 



1. For C < -2 the polynomial w(z) is negative. Thus, there are no real solutions 
ofEq. 0. 

2. For C = -2 the polynomial w(z) = -(z - l) 2 (z + l) 2 (z 2 + 2) is nonpositive, 
and it has two zeros at z = ± 1 - Setting z = ± 1 yields the singular solution 
©• 

3. For C € (-2, 0) the polynomial w(z) has 4 real roots that can be found using 
Cardano's formulae. Denoting 

, '1 . /|C|\\ /l / |C| 

a = 2 sin - arcsin — II. o = 2 cos - arccos 



3 \2 J) \3 \ 2 

c = 2 cos -arccos — , (8) 



3 \2, 

we obtain w(z) = (z 2 - a)(b - z 2 )(z 2 + c). Here < a < 1 < ft < V3 and 
c€(V3,2). 

4. For C = one obtains Schuster's solution In this case w(z) = -z 2 (z 4 -3). 
Details of the derivation of Q can be found in J5J. 

5. For C € (0,2) the polynomial w(z) has 2 real roots. Its factorisation yields 
w(z) = (z 2 + a)(z 2 + fc)(c - z 2 ), where a, ft, and c are given as before (they 
also fall in the same ranges). 

6. For C = 2 the polynomial w(z) can be easily factored, i.e., w(z) = -(z 2 - 
2)(z 2 + l) 2 . This case leads to Srivastava's solution §5§. 

7. For C > 2 the polynomial w(z) has 2 real roots. In this case the factorisation 
is even simpler. We have w(z) = (f - z 2 )(z + fz 2 + f 2 - 3), where 



f = A + l/A, A = [\(c- VC 2 -4)V , (9) 



2 

but we will not need this form of w(z). 

Out of the above 7 cases, only 3, namely (3), (5), and (7), require a separate dis- 
cussion. A closer inspection (see Sec. 2.3) shows that solution ([6]> of Goenner and 
Havas is, in fact, valid in cases (5) and (7). In case (5) a formula based on Jacobian 
elliptic functions can be also provided. 
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2.1 Solution for -2 < C < 

For C € (-2, 0) Eq. (JTJ) can be written as 

dz 1 



dt, 



^{z 2 -a){b-z 2 ){z 2 + c) 2V3 

where < kl < (here a, and c are given by Eq. Substituting x 
zl Vz 2 - a, \x\ € ( yb/ - a, oo), we obtain 



r * =-r 

J a/(z 2 - flUfe - z 2 ) (z 2 + J 



a/(z 2 - a) (i> - Z 2 ) (z 2 + cj J yj(b - a)x 2 - b a/(« + c)x 2 - c 

It is convenient to introduce a rescaled variable y = yjb - ax/ V», |v| 6(1, oo). This 
yields 

1 f dy 

where k 2 = (b - a)c/(a + c)/b. Noting that iPTOl 



/ = — 



arcdc(x, A:) = ; 1 < x < oo, 

Ji Jy^X^y 2 ^!? 

and returning to original variables, one finds a solution to Eq. ([2]) in the form 



2 V 3 V w ' V(a + c)ft 



(10) 



The new integration constant 5 corresponds to the scaling symmetry of Eq. ([2]), 
and dc is a subsidiary Jacobian elliptic function (in the standard Glaisher notation). 



Note that |dc(A\ k)\ > 1 for k e [0, 1], and thus solution ( 10) has no zeros 



dz 



2.2 Solution for < C < 2 

For C € (0, 2) Eq. leads to the integral 

7 "/ ^{z 2 + a){z 2 + b){c-z 2 ) 

where z e (— a/c, Vc). It can be solved in a similar way. The substitution y - 

V(a + c)/a(z/ Vc - z 2 ), y e R yields 

/= 1 f *! 

where again & 2 - (b—a)c/(a+c)/b. The above integral gives rise to another elliptic 
function, namely ifTOl 



arcsc(.x, Ac) = — , 



-oo < A < oo. 
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The solution can be written in the form 



acy z 



2g (ay 2 + a + c) ' 



y - sc 



(a + c)b 



ln(B£>, 



(b - a)c 
(a + c)b 



(11) 



Here again a, b, and c are given by ([8]), B is an integration constant, and sc denotes 
another subsidiary Jacobian elliptic function. Note that the expression for 6>(£) is 
regular, even for y —> ±00. The sign in (jTTJ has to be changed as the solution curve 
passes through a zero. 



2.3 Solution for C > 2 

For C > 2 and / given by Eq. ([9]), the integral 

dz 



-f 



y/(f-z 2 )(z 4 + fz 2 + f 2 -3) 



with z € (- yff, yjf) can be computed in a manner similar to that used in preceding 
sections, that is by removing the term (/ - z 2 T 1 ^ 2 from the integrand. However, it 
turns out that the substitution z = ± ^/C/(s - 1), s e [1 + C/f, 00) in the expression 



dz. 



V-z 6 + 3z 2 + C 
yields an even simpler result. One obtains 

ds 



V4s 3 - \2s - 4(C 2 - 2) 



Equation Q can can be now integrated by noticing the following identity for the 
Weierstrass elliptic function p(z;g s ,g3) with invariants g2 and #3 ifTOl : 



ds 



Returning to the original variables, we can write the solution in the form 

Vc 

y& = ± . =, (12) 

^|yp(ln(B^)/(2 V3); 12,4(C 2 - 2)) - 1 

where B is an integration constant. This result is actually valid for all C > (it 
includes solution ([TT])), and it is equivalent to solution ([6]) obtained by Goenner and 



Havas in [4 ]. The latter follows straightforwardly from the identity 

p{Xz;A~ A g 2 ,r 6 g^ = A~ 2 p(z;g 2 ,g3), 
which holds for any constant ^1^0 IfTOl . 
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Figure 1: Plot of solutions 0(£) for C = -2 (thick line) and C = -1.5, -1, -0.5. 
Only positive solutions are presented. The radial variable £ is plotted in the loga- 
rithmic scale. 



3 Discussion 

Solutions of Eq. ([2]), or equivalently Eq. (|7]), belong to three, essentially different 
types. Solutions for C € [-2, 0) are singular at £ = 0, and have no zeros. Schuster's 
integral, that is the solution for C = 0, is regular everywhere, and its sign is also 
definite. Solutions for C > are all oscillating, and the frequency of the oscillations 
grows infinitely as £ — > 0. This can be easily demonstrated for the Srivastava's 
solution Q — it has zeros at £ = exp(2kn), k e Z. 

More interestingly, all solutions for C e (-2, 0) and C > have a non-trivial, 
discrete scaling symmetry, which follows from the periodicity of elliptic functions. 
We will demonstrate this behaviour for solutions ([TO)). Applying the scaling trans- 



formation S to ( 10) one obtains 



aby 2 . 



where 



y A = dc 



(a + c)b , . . \{b-a)c 
2 \ 3 V (a + c)b 



The function dc(x, k) is periodic, with the period equal to 2K(k), where K(k) is the 
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complete elliptic integral of the first kind. Thus, for 

4V3/ 



A = exp 



V1I 



V(a + c)b 



-.K 



| (b — a)c 
(a + c)b 



m e 



(13) 



we have y& =y\, and 6(g/A)/ = Clearly, solutions ( 10) are fixed points of 



S for the values of A given by ( fT3[ ). The same formula holds for solutions ( [TT] ), and 
an analogous one can be obtained for Eq. ( fl2] ). 

Sample solutions for C < and C > are plotted in Figs.[T]and|2| respectively. 
Note that the £ axis is logarithmic on both graphs. Figure [T] shows the scale invari- 
ant solution Q — it is depicted with a thick line, and solutions given by Eq. ( 10 1 for 
C = -3/2, -1, and -1/2. The integration constant B appearing in (fT0|) is chosen in 



the way assuring that all solutions intersect at £ = 1/2. Figure|2]shows the solution 



of Srivastava (thick line), solutions given by Eq. ( |TT| ) for C = 1/2, 1, and 3/2 (thin 
solid lines), and solutions given by Eq. ( fT2] ) for C = 5/2, 3, and 7/2 (dashed lines). 
Here the constant B is set to B - 1. 

Newly obtained solutions can be used in composite stellar models on the same 
footing as Srivastava's solution [7]. Lane-Emden equations appear also in the 
description of the steady state of some field- theoretical models. In this respect, a 
good example is provided by works dealing with semilinear wave equations with 
focusing power-law nonlinearities — see e.g. lfTTl[T2l[T3l . 
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Many variants and generalisations of Lane-Emden equations appear in dif- 
ferent areas of physics (a large collection of relevant references is provided in 
[4]). These generalised versions were studied by many authors, and analytic so- 
lutions have been obtained as well. Recent literature on this subject includes also 

HUES] OH. 
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